
************************************************************************************************************************************************
************************************************************************************************************************************************
*** Paper "Drivers of Political Participation: The Role of Partisanship, Identity and Incentives in Mobilizing Zambian Citizens"   
*** Replication Dofile for Submission to Comparative Political Studies
*** Authors: Prisca Jöst, Matt Krönke, Sarah Lockwood and Ellen Lust
************************************************************************************************************************************************
************************************************************************************************************************************************


************************************************************************************************************************************
******* Robustness Check: New Analysis using Data after merging with information on party strongholds (Table D6 in Appendix) ******* 
************************************************************************************************************************************

*** We merged experimental data with data on party strongholds and change in partistrongholds between 2016-2021
*** Not merged: We lost 334 respondents in constituencies for which stronghold data was not available 
*** Sub-sample analysis: note dropping decisions

clear all

set maxvar 30000

use "Const_Data_For_Replication.dta"

generate participation = .
replace participation = 1 if q1y == 4
replace participation = 2 if q1y == 3
replace participation = 3 if q1y == 2
replace participation = 4 if q1y == 1

* rescale variable to |0,1]
generate participation2 = (participation - 1) / (4-1)

* generate CoGender = .
gen CoGender = .
replace CoGender = 1 if demo_q4 == 2 & Gender == 2 | demo_q4 == 1 & Gender == 1
replace CoGender = 0 if demo_q4 == 2 & Gender == 1 | demo_q4 == 1 & Gender == 2

tab CoGender


drop if C1 == "Don't Know/Refuse to Answer"

** For Patriotic Front Only **

gen PF = 0
replace PF = 1 if Stronghold1 == "PF" & Stronghold2 == "PF"

keep if PF==1

eststo MainStr1 : reg participation2 i.Activity i.Authority ib2.CoEthnic ib2.Partisanship ib2.Origin ib2.Residence i.CoGender ib2.SocialBenefit ib2.Payment i.SanctioningLeader i.SanctioningCommunity

esttab MainStr1, se star(+ 0.10 * 0.05 ** 0.01 *** 0.001)
*esttab using BasicStr1.doc, se star(+ 0.10 * 0.05 ** 0.01 *** 0.001)

coefplot, xline(0) keep(*CoEthnic *Partisanship *CoGender *SocialBenefit *Payment *Origin *Residence) name(MainStr1)	

** For United Party for National Development Only **

gen UPND = 0
replace UPND = 1 if Stronghold1 == "UPND" & Stronghold2 == "UPND"

keep if UPND==1
eststo MainStr2 : reg participation2 i.Activity i.Authority ib2.CoEthnic ib2.Partisanship ib2.Origin ib2.Residence i.CoGender ib2.SocialBenefit ib2.Payment i.SanctioningLeader i.SanctioningCommunity

esttab MainStr2, se star(+ 0.10 * 0.05 ** 0.01 *** 0.001)
*esttab using BasicStr2.doc, se star(+ 0.10 * 0.05 ** 0.01 *** 0.001)

coefplot, xline(0) keep(*CoEthnic *Partisanship *CoGender *SocialBenefit *Payment *Origin *Residence) name(MainStr2)


** For Party Stronghold Changed 2016-2021 Only **

keep if Change==1
eststo MainStr3 : reg participation2 i.Activity i.Authority ib2.CoEthnic ib2.Partisanship ib2.Origin ib2.Residence i.CoGender ib2.SocialBenefit ib2.Payment i.SanctioningLeader i.SanctioningCommunity

esttab MainStr3, se star(+ 0.10 * 0.05 ** 0.01 *** 0.001)
*esttab using BasicStr3.doc, se star(+ 0.10 * 0.05 ** 0.01 *** 0.001)

coefplot, xline(0) keep(*CoEthnic *Partisanship *CoGender *SocialBenefit *Payment *Origin *Residence) name(MainStr3)
	
